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The contribution of the Warm-Hot Intergalactic Medium 
to the CMB anisotropics via the Sunyaev-Zeldovich effect 
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ABSTRACT 

Cosmological hydrodynamical simulations predict that a large fraction of all baryons 
reside within mildly non-linear structures with temperatures in the range 10 5 — 10 7 K. 
As the gas is highly ionized, it could be detected by the temperature anisotropies 
generated on the Cosmic Microwave Background radiation. We refine our previous es- 
timates of the thermal Sunyaev-Zeldovich effect by introducing a non-polytropic equa- 
tion of state to model the temperature distribution of the shock heated gas derived 
from temperature-density phase diagrams of different hydrodynamical simulations. 
Depending on the specific model, the Comptonization parameter varies in the range 
10~ 7 s$2/ c <2xlCr 6 , compatible with the FIRAS upper limit. This amplitude is in 
agreement with a simple toy model constructed to estimate the average effect induced 
by filaments of ionized gas. Using the log-normal probability density function we calcu- 
late the correlation function and the power spectrum of the temperature anisotropics 
generated by the WHIM filaments. For a wide range of the parameter space, the max- 
imum amplitude of the radiation power spectrum is {£ + l)lCt/2ir = 0.7 — 70(fiK) 2 
at £ « 200 — 500. This amplitude scales with baryon density, Hubble constant and the 
amplitude of the matter power spectrum og as [(£+ l)£Ct] max /2ir oc of' 6 '(ilbh) 2 . Since 
the thermal Sunyaev-Zeldovich effect has a specific frequency dependence, we analyze 
the possibility of detecting this component with the forthcoming Planck data. 
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1 INTRODUCTION 

The distribution of baryons in the local Universe is one of 
the main problems of modern Cosmology. The highly ion- 
ized intergalactic gas, that evolved from the initial density 
perturbations into a complex network of mild ly non-linear 
structures in the r edshift interval 2 < z < 6 (|Raucb|[l998l ; 
IStocke et al. 1 120041), could contain most of t he baryons in 
the u niverse ( Rauch et al.lll997l ; ISchave]|200ll ; lRichter et al.l 
2006). At redshifts z > 2 most baryons are found in the Lya 
systems detected through absorption lines in the spectra 
of distant quasars. With cosmic evolution, the baryon frac- 
tion in these structures decreases as more matter is concen- 
trated within compact virialized objects. At low redshifts, 
the Lya systems are filaments wit h low HI column d ensities 
containing ~ 30% of all baryons (|Stocke et al.l [20041 ). while 
the material in stars and galaxies contains about 10% of all 
baryons. An extra 5% coul d be in the form of Circumgalactic 
Medium around galaxies l|Gupta et al.|[2012l ). More impor- 
tantly, at z ~ about 50% of all cosmic baryons have not yet 
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been identified jFukugita fc Peebles 2004; Sh ull et alJlioTlh . 
Hydrodynamical simulations predicted that baryons could 
be in the form of shock-heated intergalactic gas in mildly- 
nonlinear structures with temperatures 10 5 — 10 7 K, called 
Warm-Hot Intergalactic Medium (WHIM). The baryon frac- 
tion in this medi um could be as larg e as 40%-50% in 
the local Universe JCen fc Ostriker 199S ; Dave et al1 ll999 



200ll;lDanforth fc ShulfcOOSMProchaska fc Tumlinsonlbood ; 
Smith et al.ll201l| j containing the bulk, if not all, the uniden- 
tified baryons. 

The observational effort has concentrated in searching 
for the WHIM X-ray s ignature both in emission and in ab- 
sorption. ISoltanl (|2006t ) looked for the extended soft X-ray 
emission around field galaxies but his task was complicated 
by the need to subtract all systematic effects that could 
mimic the diffuse signal. More recently, the effort has been 
concentrated in searching for absorption lines due to highly 
ionized heavy elements fr om the far-ultraviolet to the soft X- 
ray fsee lShull et alj 2011, for a revi ew) with partial success. 
For example. iDietrich et al.l (|2012l 'l have recently reported 
the detection of a large-scale filament between two clusters 
of galaxies. 
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While absorption lines will identify individ- 
ual systems, complementary techniques are needed 
to study the overall properties of the WHIM . To 
this purpose, lAtrio-Barandela fc Miicketl ((2006) and 
lAtrio-Barandela et alj l|2008h suggested that the WHIM 
distribution would generate temperature anisotropics on 
the Cosmic Microwave Background (CMB) due to the 
therma l (TSZ) and kinematic (KSZ) S unyaev-Zeldovich 
effect (|Sunvaev fc Zeldovichl 1 19721 , |l98Ch . opening a new 
observati onal window to study the g as distribution at low 
redshifts (|Genova-Santos et al'j|2009h . lHallman et alj (|2007T ) 
found that about a third of the SZ flux would be generated 
by unbound gas, but their estimated amplitude was lower 
than the prediction of our model. Our computation as- 
sumed that the weakly non-linear matter distribution was 
describ ed by a log-normal p r obability distribut i on fun ct ion 
PDF (|Coles fc Jones! Il99ll ; IChoudhurv et al.1 l200lh and 
the gas followed a polytropic equation of state. While the 
IGM at redshifts z > 2 behaves like a polytrope, at low 
redshifts the polytropic equation of state is not longer valid 
since the gas is heated mostly by shocks to temperatures 
about 10 5 K < T < 10 7 K at density contrasts 8 < 100. 
In the regions where the gas is shock-heated, small-scale 
density perturbatio ns can be neglected due to the high 
pressure gradients l|Klar fc Miicketl |2010| V simplifying the 
treatment. In this paper we shall improve our previous 
calculations by using equations of state derived from 
hydrodynamical simulations to provide a more physical 
description of the WHIM. Briefly, in section [2] we estimate 
the order of magnitude contribution of a population of 
WHIM filaments to the Comptonization parameter based 
on a crude geometrical model; in section [3] we describe our 
log-normal model and we detail the differences with our 
previous treatments; in section [4] we compute the CMB 
Comptonization parameter due to the WHIM and in section 
[5] the correlation function and angular power spectrum; in 
section [6] we present our results and, finally, in section [7] we 
summarize our conclusions. 



2 ESTIMATE OF THE MEAN 

COMPTONIZATION PARAMETER DUE TO 
THE WHIM 

The temperature anisotropies due to the inverse Compton 
scattering of CMB photons by the free electrons have only 
been measured for clusters of galaxies. There are two contri- 
butions: the thermal component (TSZ) due to the motion of 
the electrons in the potential wells and the kinematic com- 
ponent (KSZ) due to the motion of the cluster as a whole. 
The induced temperature anisotropies are 



JO /TSZ \ JO /KSZ 



Vcl_ 

c 



(1) 



where G(x) = x coth(i/2) - 4, x = hv/k B T is the CMB fre- 
quency in dimensionless units, Vol is the velocity of the clus- 
ter projected along the line of sight (los) and c the speed of 
light. The Comptonization parameter and the optical depth 
r are defined as 



V = yo / T e n e dl, 



n e dl , 



(2) 




Figure 1. Geometrical model of a simple filament with constant 
density and temperature. The filament has a size L, width d, is 
located at redshift z at a comoving distance -D(z), subtends an 
angle a from the observer and is randomly oriented in space, 
forming an angle 8 with respect to the line of sight. 



where T e and n e are the electron temperature and electron 
number density, dl is the proper distance along the los and 
yo = crrkB /m e (? , with <jt, ks, m e Thomson cross section, 
Boltzmann constant and electron mass, respectively. While 
the KSZ effect has the same frequency dependence as the 
intrinsic CMB anisotropies, the change with frequency of 
TSZ effect, G(x), is different from that of any other known 
foreground. It is negative in the Rayleigh- Jeans and positive 
in the Wien regions, being null close to v — 217GHz and can 
be more easily separated from other contributions than KSZ 
anisotropies. Since they are more easily detectable, here we 
shall discuss only the TSZ contribution of the WHIM. 

An order of magnitude estimate of the WHIM TSZ ef- 
fect can be obtained by considering the effect on the CMB 
radiation of a single filament with constant electron density 
n e and constant temperature T e . The geometrical configura- 
tion of a filament of size L and thickness d is shown in Fig- 
ure [1] For a filament randomly oriented u = (8, <j>) in space, 
let 6 be the angle between the filament and the los. The 
Comptonization parameter induced in the radiation cross- 
ing this filament is 



V = yo- 



, T e d 



cos 9 



(3) 



where d/ cos8 is the projected thickness of the filament. 
The Comptonization parameter averaged over all possible 
azimuths d> is 



(y)< 



yo- 



He Ted d 

cos 2nL' 



(4) 



If the size of filaments at redshift z is small compared to its 
angular diameter distance, L <g Da (z), then the projected 
comoving length L' can be written as L' = LcosO ~ aD(z) 
where D(z) is the co-moving distance and a is the angle 
subtended by the filament on the sky. Then, 



{y)<t> = yo 



n e T e Ld 2 
2%a 2 D 2 (z) 



(5) 



The probability that a given los crosses the filament is 
L/(nD(z)); if d <g L, then the average over all possible 
values of a is 



(y) = yo 



Tie Te L 

2tt 2 Djz) 



d. 



(6) 
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In their ana lysis of the properties of filaments, 
iKlar fc Mucketl (2012h obtained d = eL, with e G [0.01-0.1]. 

The mean distortion along the los can be obtained by 
adding the effect of all filaments of different masses at dif- 
ferent redshifts: 



(y) 



_L_ 

2^ 



dnf 
dM 



dMy n e T t 



eL 2 (M) dV(z) 



dz 



dz , 



where rif is the number of filaments per unit of volume. To 
proceed further we assume that filaments connect clusters 
and groups of galaxies. Then, their length is approximately 
the distance between clusters. If n c ; is the number of clusters 
of a given mass per unit of volume, then the mean separation 
of the filaments within a volume V is I ~ (V/nd) 1 ^ 3 and its 
mean length L ~ n cl . Thus, the mean Comptonization 
parameter along the los is 



(y) 



with 



dnf 
dM 



dM y n e T e e 



D(z) dD{z) 



n c i(M)' 2 / 3 dz 



dz , (8) 



dD(z) = aff _1 [fiA + fim(l + zf + Qfe(l + z) 



2,1/2 



dz . 



(9) 



and Qfc = (1 — fl m — Qa)- Equation © can be used to ob- 
tain an order of magnitude of the distortion induced by the 
electrons on the WHIM filaments. If we take T e « 10 7 K, 
e = 0.05 and compute the number of filaments according to 
Sheth & Tormen (2002) then y c = 10" 6 , compatible with 
the FIRAS upper limit of y c = 1.5 x 10" 5 . The CMB distor- 
tion due to a network of filaments of ionized gas is compat- 
ible with observations. The integration is performed from 
Mi = 10 13 M Q to M f = 10 15 M Q . Our results are insensitive 
to the exact value of Mf whose number density is exponen- 
tially suppressed but depend on Mi which we took to be 
the mass of rich groups of galaxies, limit where the number 
of objects is not longer well described by Sheth fc Tormen 
(2002). Using that formalism IShimon et al.l (|2012T ) have cal- 
culated the radiation power spectrum. In the next section 
we will give a more refined prediction based on a log-normal 
model of the baryon distribution. 



3 THE LOG-NORMAL BARYON DENSITY 
DISTRIBUTION MODEL 

In perturbation theory, the density contrast is defined as 

P(x)-(P> 



f(x) 



(10) 



where p(x) is the density at any given point and (p) the 
mean density at any given redshift. In Fourier space and in 
the linear regime, all modes evolve at the same rate and 



<5(k,2) = D + (z)5(k,0) . 



(11) 



where D+(z) = D+(z, Ha, fim) is the linear growth factor 
normalized to D+(0) = 1. If the initial density field is a 
Gaussian random field, then the linearly extrapolated power 
spectrum P(k) is defined by 

(5(k, 0)«5(k', 0)) = (27r) 3 P(k)5(k - k') . (12) 

In eq.[T2]the power spectrum P(k) is that of the background 
ACDM model, computed using linear theory. We took the 
spectral index at large scales to be ng = 1. We did not 



consider corrections due to the non-linear evolution of the 
matter density field on scales about and below 8/i _1 Mpc. 

In the linear regime, baryons follow the Dark Matter 
(DM) distribution, but this is not longer true once struc- 
tures become non- linear . In regions wher e large -scale pertur- 
bations develop shocks, IKlar fc Miicketl l|2010f ) showed that 
the gas perturbations on small scales were suppressed by the 
enhanced pressure and would be nearly erased. Then, the 
distribution of the baryon WHIM filaments that are form- 
ing and evolving within the large-scale density perturbations 
would be different from that of the DM. The linear density 
contrast of baryons in the IGM can be obtained from that of 
the DM by smoothing over scales belo w Lp, the la rgest scale 
erased by shock- heating. According to iFang et all l| 19931 ) 



5b(K z) 



5dm (k, z) 
1 + L%k 2 



(13) 



The comoving length scale Lq can be determined by im- 
posing that it is the smallest possible length scale at which 
the linear peculiar velocity v p equals the sound speed v s (z) 
of the baryon fluid. The sound speed is determined by the 
mean temperature Tigm(z) of the IGM at any given z: 
v s = y // 2kBTiQM{z)/m p . The condition v p ^ v s gives 



Lo(z) 



2tt(1 + zf VsH- 1 



(14) 



(n A + n m (i + z)3y/ 2 D + (z)5 ' 

In eq. (|14|l . the sound speed is determined by the IGM 
temperature Tigm at mean density. Tigm ~ 10 4 if and it 
is mainly determined by the evolution of the UV back- 
ground. For that temperature, at redshift z — we have 
Lo ~ \.7 h~ x Mpc. This value of cut-off scale agrees with the 
results of lKlar fc Miicketl (120121) on the formation and evo- 
lution of WHIM filaments~ TTittlev fc Meiksinl (120071 1 have 
studied different ionization models and discussed the pos- 
sible range of variation of the mean IGM temperature. In 
our calculations, we assume that Tigm varies as given in 
iTheuns et all (|2002h . Below we shall see that at redshifts 
z > 3 the WHIM does not generate significant anisotropies. 
At redshifts 2^3, the mean IGM temperature variation 
with redshift is small and can be roughly approximated as 
log 10 (TiG M /10 3 K) w (A + 0.1(1 + 0)), where A varies in the 
range A = 0.5 - 0.9, corresponding to Tigm = 10 3 6 - 10 4 K. 

The Fourier transform of eq. (|13l) gives the density con- 
trast in real co-moving space, 5(x, z). If the density field 
grows non-linearly but the velocity field still remains in the 
linear regime, then the initial Gaussian den sity field is de- 
scrib ed by a log-normal density distribution (I Coles fc Jonesl 
1991). This log-normal model has been shown to de- 
scribe the matter distribu tion in hydrodynamical simula- 
tions (|Bi fc Davidsen 1997[) and to repro duce the observa- 
tions of the IGM (jChoudhurv et al.ll200ll ). For a log-normal 
density field, the baryon number density is given by 



tib(x, z) = m(z)e 



} W/2 e i fl (x,«) 



(15) 



where 8b is the linear baryon density field and A% = (6%) 
its variance, with 



A b (z) = {5 2 b(x,z))=D 2 + (z) 



P DM (k) d 3 k 



(16) 



[1 + fc 2 Lg]2 (2tt) 3 

If 5b is Gaussian distributed, then the spatial mean 
{iib(x,z)) equals the cosmic baryonic background density 
no(z). Hereafter we shall denote by £ = w_b(x, z)/no(z) the 
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non-linear baryon density in units of the baryon mean den- 
sity. 

In lAtrio-Barandela fc Miicket] J2006I) 



lAtrio-Barandela et alj ()2008l ') and iGenova-Santos et al.l 
(|2009h we assumed that the gas in the WHIM follows 
a polytropic equation of state. High resolution hydro- 
simulations have shown that at the physical conditions 
of the WHIM, the shock heated gas at temperatures 
ltfK < T e < W 7 K and den s ity co ntrast s 1 < £ 100 
is not polytropic. iKang et al.l (|2005l ) and ICen fc Ostrikerl 
(2006) provide phase diagrams of electron temperature T e 
versus baryon density ub in the range of interest. Since 
the IGM is highly ionized due photo-ionization by the UV 
background, then n e ~ lis. This is even more accurate for 
the WHIM as collisional ionization at high temperatures 
also contributes to maintain n e = ns- We have used these 
results to construct fits to the baryon temperature- density 
distribution, T e = T e (n e ). We used the fitting formula 



log'i 



10 8 K 



log 10 (4 + £ Q+1 /«) 



(17) 



This expression fits the results of IKang et al l (|2005h when 
a = 1 -4. 

The phase space derived from simulations are not sim- 
ple linear relations of the type T e — T e (n e ); phase diagrams 
have a large scatter. This theoretical uncertainty is impor- 
tant since the Comptonization parameter and temperature 
anisotropies depend on how steeply the temperature rises 
with increasing density contrast. We model this uncertainty 
by taking a = 1 — 4 in eq. (|17[) , a range wide enough 
to reproduce sufficiently well the phase diagrams and the 
variation therein. We have also considered the results of 
ICen fc Ostrikerl l|2006l ); after correcting a possible typo their 
expression reads log 10 (T e /10 8 K) = -2.5/ log 10 (4 + £) a9 . 
The phase diagrams obtained from c osmological simu l ations 
agree qualitatively with the results of lKlar fc Miicketl (|2012l ) 
on the effect of shocks on WHIM filaments. In particular, 
as time evolves shocks propagate to lower density regions, 
steepening the T e — n e relation as shown in the phase dia- 
grams of numerical simulations. 

For illustration, in Fig.[2]we plot a series of different fits. 
The li nes with labels T3 and T4 correspond to IKang et al.l 
(2005) m odel with a = 1 .5 an d 4, respectively, T2 corre- 
sponds to lCen fc Ostrikerl (|2006l ) with the corrected formula 
given above. For comparison, Ti corresponds to a simple 
power law with T e oc n^ 2 . From Fig. [2] we can immediately 
conclude that since filaments of any density in the T3, T4 
models are hotter than in the T\,Ti models, temperature 
anisotropies and distortions would be larger in the former 
than in the latter cases. 



4 THE COMPTONIZATION PARAMETER 

The log-normal formalism introduced in the previous section 
permits to compute the statistical properties of the TSZ 
generated by the baryon population of WHIM filaments. 
The average Comptonization par ameter y av along a los is 
ijAtrio-Barandela fc Miicketll200r3 ) 




2/av = Vo / (n e (z,x)T e (z,x))dl , 



(18) 



100 



Figure 2. Various fit functions T e = T e (n e ) to phase diagrams 
derived from hydro-numerical simulations. Curves T3, T4 corre- 
spond to the results of Kang et al. (2004) with a = 1.5, 4, re- 
spectively; T2 corresponds to Cen &; Ostriker (2006); finally, Ti 
represents the power law T oc n 3 / 2 . 



where the electron and temperature distribution are given 
by eqs. (|15[) and (|17|l . respectively. 

The average in eq. (|18p must be limited to those scales 
(£1,62) where the non-linear evolution is well described by 
the log-normal probability density function. The choice of £1 
is somewhat arbitrary. Since only gas with density contrasts 
£ > 2 will undergo shock-heating, this could be taken as a 
lo wer lim it. As a criteria, we take £1 = £median+ff, where a — 
yj uar(£) (war(£) = exp(Aa) — 1) is the standard deviation of 
the log-normal distribution . To include only the scales that 
undergo shocks we also require that £1 > 2 at all redshifts. 
As upper limit we took £2 ^ 100, but the results were weakly 
dependent on this value. 

Even though this criteria is somewhat arbitrary, the 
mass fraction on the interval (£i,£2) is very close to that of 
the undetected baryons: 40-50%. At any redshift, the mass 
fraction of the gas contained in filaments with overdensities 
in the range [£1 , £2] depends on the IGM temperature. If £ is 
the mean of the distribution then M(£i,£ 2 ) = £ _1 £ F (0 
is the fraction of mass in the integration range. At z = 0, 
the fraction of baryons in the WHIM is « 48% for A = 0.9 
and « 43% for A = 0.5. One would expect the mass frac- 
tion to be largest for Tigm = 10 3 ' 6 K, since smaller scales 
would survive shock heating than for Tigm = 10 4 K. The 
difference comes from the lower limit £1. For A — 0.5, the 
rms of the log-normal distribution a is larger and so is £1 ; 
the integration is carried over a smaller range of densities 
and low density regions that are included for A = 0.9 are 
excluded for A = 0.5. However, even if in this aspect the 
model is not fully satisfactory, the low dense regions do not 
contribute significantly to spectral distortions or tempera- 
ture anisotropies and the exact value of the lower limit £1, 
while important for the baryon fraction on the WHIM, has 
little effect on the final results. 



5 THE ANGULAR POWER SPECTRUM 

The correlation function of the CMB temperature 
anisotropies due to the WHIM TSZ contribution along two 
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different los 1 and 1' separated b y an angle 8, is given by 
l|Atrio-Barandela fc Mucketll2006h 

C(9) = y 2 J^ J\s{^z)S^',z l ))f z ^ ; dzdz' , (19) 

where S(x.,z) = n e (x, z)T(n e (x, z)). The density average of 
eq. (| ltjp is carried out over the same density interval than in 
eq. flg). 

Since the equation of state of the gas is not longer poly- 
tropic, the correlation of the electron density at different 
positions can not be factored out. The average in eq. (|19[) 
differs from our previous treatment since now it needs to be 
computed using the bi-variate log-normal PDF. Remember- 
ing that £ = iib(x, z)/no(z), then 



(S(x,z)S(x',z'))= I I tT(m'T(ti')F(i,t')dSdt f , (20) 
where the bi-variate log-normal PDF F(^,^') is given by 



exp 



2(l-r^) 



7T X 



(logg-M) 2 9 „ (log|-M)(log C'-m') , (logg'-M')^ 

in ^"c a _ a > 1 



(21) 



In this expression, As is given by eq. ()16[) . = — A fl /2 is 
the mean of the log-normal PDF and r c is the correlation 
coefficient r c = (log(£) log(^'))/(As A' B ) that in this partic- 
ular case is given by 

_D+(z)D+(z') f P DM (k)j () (k\yL-yL'\)k 2 dk 



27r 2 A B A' 



[1 + L {z') 2 k 2 ][l + L (z) 2 k 2 



(22) 



In this expression, jo denotes the spherical Bessel function 
of zero order and |x — x'| the proper distance between two 
patches at positions x(z) and x(z') separated an angular 
distance 6. In the flat-sky approximation 



|x - x'| a ^l ± (6,z)z + [r(z)-r(zi)]2, 



(23) 



with l±(0, z) the transverse distance of two points located at 
the same redshift. Therefore, the angular dependence of the 
correlation function in eq. (|19[1 enters only through eq. (|23[) 
in eq. fl22j. 

The power spectrum can be obtained by the Fourier 
transform of the correlation function (eq. I19|l 



Co = 2tt 



+ i 



C(6)P e (cos(6))d(cos(6)) 



(24) 




where Pi denotes the Legendre polynomial of multipoles 



redshift z 

Figure 3. Comptonization parameter 10 6 y as a function of red- 
shift. From top to bottom solid and dashed lines corresponds to 
model T4 and A = 0.5, 0.9, dotted and dot-dashed lines to T3 
and A = 0.5, 0.9, respectively. Finally the double dot-dashed line 
corresponds to T2 with A = 0.5. 



characterized by a parameter A w 0.5 — 0.9 (see sec. [3} and 
for the latter we shall consider three different relations: A 
polytropic equation of state T e oc rfj 2 , the Cen & Ostriker 
(2006) model and Kang et al (2004) model given by eq. iflS) 
with a — 1.5, 4, denoted by by Ti, T2, T3, T4, respectively. 



6.1 The mean Comptonization parameter 

In Fig. [3] we plot the contribution to the average Comp- 
tonization parameter y av from z — up to the given redshift 
for different models. In all cases, most of the contribution 
comes from z ^ 1. In retrospect, this justifies restricting 
in sec. [3] the IGM temperature parameter to the interval 
A = 0.5 — 0.9, valid for the interval z ^ 3, since baryons at 
higher redshifts do not contribute to the TSZ effect. 

In the figure, the amplitude of the CMB distortion de- 
pends on the equation of state and the temperature of the 
IGM. When T IG m = 10 3 6 K smaller (and therefore more) 
scales contribute to the shock-heated WHIM and Lo is 
smaller. Those small scales contribute to the variance of the 
baryon field at all redshifts (see eq. I16|l . These regions are 
of high density and at high temperature and the total ef- 
fect increases, compensating for the baryon fraction being 
smaller (see section |4}. 



6 RESULTS 

To compute the temperature anisotropies generated by the 
WHIM we take as a background cosmological model the 
concordance ACDM model with densities dark energy, dark 
matter and baryon fractions J1a = 0.75, fl m = 0.25, fij, = 
0.043 and Hubble constant h = 0.73 and the amplitude of 
the matter density perturbations at 8ft _1 Mpc, ug = 0.8. 
For simplicity, we fix the frequency dependence of the TSZ 
effect to unity, G(v) = 1. Our results depend also on physical 
parameters such as Lo and the shape of the equation of state. 
The former is determined by the mean IGM temperature, 



6.2 The correlation function and the power 
spectrum 

In Fig. [4] we show the power spectrum of the radiation tem- 
perature anisotropies generated by the WHIM for different 
temperature models, obtained by numerical integration of 
eqs. (|19I24|I . Lines follow the same convention than in Fig. [2] 
In Fig. Owe show the dependence with equation of state and 
Tigm- The overall shape of the power spectra are very sim- 
ple: the functional form has a single maximum. For the range 
of variation of the model parameters, the maximum occurs 
in the range £ = 200 — 500 and the maximum amplitude 
is £(£ + l)C e /2n ~ 0.07 - 70{fj,K) 2 . The different equations 
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Figure 4. Power spectrum of the radiation temperature 
anisotropies generated by the WHIM. Lines correspond to the 
models given in Fig[2] with A = 0.5. 



Figure 6. Contribution to the power spectrum of different red- 
shift intervals, for the T4 model with A = 0.5. 
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Figure 5. Effect of the cut-off length in the radiation power 
spectra. Curves Mi, M2 correspond to A = 0.5, A = 0.9, respec- 
tively. From top to bottom, the equation of state corresponds to 
the T 4 , T 3 and T 2 models. 

of state Ti — T4 give power spectra that differ by three or- 
ders of magnitude in amplitude. This wide range reflects our 
theoretical uncertainty on the equation of state of the shock- 
heated gas. Using a polytropic equation of state gave even a 
larger range, and a better determination of the TSZ power 
spectrum can only come through better understanding of 
the physics of the WHIM, that is, through cosmological hy- 
drosimulations with larger dynamical range. 

The effect of the IGM temperature changes both the 
amplitude and the location of the maximum. This was to 
be expected since for smaller IGM temperature the cut-off 
length Lo is smaller and the contribution of small scales 
is more important, increasing the amplitude reducing the 
average angular size of the anisotropies. For A — 0.5 the 
power spectrum reaches a maximum at £ ^ 450 with an 
amplitude of g{l + l)CV, max /27r w 70(fxK) 2 . For A = 0.9, 
the maximum is at I ~ 300 and its amplitude is smaller by 
a factor ~ 3. 

In Fig. [6] we show the contribution to the power spec- 
trum due to the WHIM at different redshift intervals. Sim- 



ilarly to what occurs with the Comptonization parameter, 
the dominant contribution comes 0.3 < 2 < 0.5, and the 
contribution for structures with z > 1 is negligible, in agree- 
ment with the resu lts of numerical simulations reported by 
ISmith et al.1 (|201lh . A gain, this justifies the use of a simple 
parametrization of the IGM temperature that is only valid 
for redshifts z ^ 3. 

The radiation power spectrum of the WHIM induced 
temperature anisotropies also depends on the cosmological 
parameters of the concordance model. As the amplitude of 
the TSZ effect depends on the number of electrons projected 
along the los, the amplitude of the effect scales as Qsh. 
Other parameters like Q m or change the induced tem- 
perature anisotropies because they change the shape matter 
power spectrum in eq . (|16[) . But the variance of the linear 
baryon density field A| is dominated by the amplitude of 
the matter power spectrum at small scales, and the largest 
effect is produced by variations in the amplitude at those 
scales. Then, the largest variations are produced by changes 
in ag. These are shown in Fig. [7] A linear regression fit to 
the power spectra of different cosmological parameters gives 
the following scaling relation at the maximum 

[(I + l)lCi] m ^/2n oc al\VL b h) 2 . (25) 

This scaling relations differ from our previous results in two 
very important aspects: the power spectrum does not de- 
pend so strongly on the model parameters and the TSZ 
contribution of the WHIM is always smaller than that of 
clusters. However, like in Atrio-Barandela & Miicket (2006), 
the bulk of the contribution occurs well within the angu- 
lar range measured by WMAP and Planck. This opens the 
possibility of looking for the WHIM contribution in the ra- 
diation power spectrum as in Genova-Santos et al (2009). A 
more suitable data set is the forthcoming Planck data. The 
high resolution and wide frequency coverage of this instru- 
ment makes it feasible to trace the WHIM not only by its 
power spectrum, but also by the changing amplitude with 
frequency. In particular, the 217GHz channel, that corre- 
sponds to the TSZ null, will help to separate the signal from 
foreground contributions and other larger systematics. 
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Figure 7. Radiation power spectrum for the T4 model with A 
0.5 and different values of erg. 



7 CONCLUSIONS 

We have shown that the gas within non-linear large-scale fil- 
aments, with density contrasts 2 sC £ ^ 100, shock-heated up 
to temperatures 10 5 K < T < 10 7 K will produced measur- 
able temperature anisotropies via the TSZ effect. We have 
assumed that the gas follows the DM distribution that is well 
described by a log-normal PDF in the density range of inter- 
est, in agreement with N-body simulations. The main differ- 
ence with Atrio-Barandela (2006) and Atrio-Barandela et al 
(2008) is the adoption of an equation of state T e = T e (n e ), 
appropriate to describe the shock-heated gas in the WHIM. 
The equation of state was derived from phase-diagrams of 
cosmological hydrodynamical simulations. This complicated 
the statistical treatment, requiring the statistical averages 
over densities to be carried out using the bi-variate log- 
normal PDF. Our final results depend on the cosmological 
models, the temperature of the IGM and the functional form 
of the equation of state. We showed that for temperatures 
compatible with observations and for functional relations 
derived from numerical simulations, the amplitude of the ra- 
diation power spectrum is £{£ + l)C e /2n = 0.07 - 70(fj,K) 2 , 
being its maximum amplitude at £ ~ 200 — 500 and with 
the bulk of the contribution coming from redshift z ^ 1.0. 
The large scatter in the final values of the Comptoniza- 
tion parameter and in the amplitude of the power spec- 
tra reflects the theoretical uncertainty associated with the 
functional shape of the temperature-density relation of the 
shock-heated gas. 

The WHIM contribution peaks at a much larger angu- 
lar scale than £ ~ 3000, that of clusters of galaxies, mak- 
ing it accessible to the obser vations of WMAP and Planck . 
Recently we have analyzed (jSuarez- Velasquez et al.l 120121 ') 
the observational prospects of the WHIM TSZ by cross- 
correlating CMB data with tracers of the gas distribution. 
We concluded that for the current model parameters, the 
density field reconstructed from the 2MASS galaxy catalog 
and Planck data could provide the first evidence of the large 
scale distribution of the WHIM. 
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